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Abstract. The Lax pseudo-differential operator plays a key role in studying the general 
set of KP equations, although it is normally treated in a formal way, without worrying 
about a complete characterization of its mathematical properties. The aim of the present 
paper is therefore to investigate the ellipticity condition. For this purpose, after a careful 
evaluation of the kernel with the associated symbol, the majorization ensuring ellipticity is 
studied in detail. This leads to non-trivial restrictions on the admissible set of potentials in 
the Lax operator. When their time evolution is also considered, the ellipticity conditions 
turn out to involve derivatives of the logarithm of the r-function. 



PACS numbers: 05.45.Y 



1 



1. Introduction 



Several important developments in modern mathematical physics are due to the investiga- 
tion of pseudo-differential operators on R m and on general Riemannian manifolds [1] . For 
our purposes, it is sufficient to recall the following basic properties. 

(i) A linear partial differential operator P of order d can be written in the form 

p=j2 a «w D * ( L1 ) 



|a|<d 



where (here i = V - 1) 



rn 



\a\=J2 a k (1.2) 



k=i 



^ s ^id n -&°" (L3) 



and a a is a C°° function on R m for all a. The associated symbol is, by definition, 

\a.\<d 

i.e. it is obtained by replacing the differential operator D% by the monomial £ a . The pair 
(x, £) may be viewed as defining a point of the cotangent bundle of R m , and the action of 
P on the elements of the Schwarz space S of smooth complex-valued functions on R m of 
rapid decrease is given by 

Pf(x) = J e^-y^p(x,Of(y)dydt (1.5) 

where the = dy\...dy m and = <i£i...<i£ m orders of integration cannot be interchanged, 
since the integral is not absolutely convergent. 

(ii) Pseudo-differential operators are instead a more general class of operators whose symbol 
need not be a polynomial but has suitable regularity properties. More precisely, let S d be 
the set of all symbols p{x, £) such that [1] 



(1) p is smooth in with compact x support. 

(2) For all (cc, /?), there exist constants C aj/ g for which 

<C a ^(l + |e|) d - |/31 (1.6) 

for some real (not necessarily positive) value of d, where |/3| = X)feLi ^ ( see (1-2))- The 
associated pseudo-differential operator, defined on the Schwarz space and taking values in 
the set of smooth functions on R m with compact support: 

P:<S^C c °°(R m ) 

is defined in a way formally analogous to Eq. (1.5). 

(iii) Let now U be an open subset with compact closure in R m , and consider an open 
subset Ui whose closure U\ is properly included into U: U\ C U. If p is a symbol of order 
d on U, it is said to be elliptic on U\ if there exists an open set U2 which contains U 1 and 
positive constants Ci so that 

\p(x,0\- 1 <C 1 (l + \t\)- d (1.7) 

for |£| > C and x G U2, where 

(i-s) 

*!=i 

The corresponding operator P is then elliptic. 

From a mathematical point of view, pseudo-differential operators occur in many prob- 
lems in global analysis [1, 2], and recent developments deal with the functional calculus 
of pseudo-differential boundary- value problems [3]. From a physical point of view, such 
a formalism is important in quantum gravity and quantum field theory [4-6]. In partic- 
ular, we are here interested in an interdisciplinary field, i.e. a rigorous approach to the 



D«D^p(x,0 



ICI = y/g ab (x)Ub 



\ 
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Kadomtsev-Petviashvili (hereafter KP) equations. Recall that the KP equation can be 
written in the form [7, 8] 



d ( ' du . „ du d s u\ „ 2° 2u 



dx \dt ^^dx^ dx 3 J ^ a dy 2 ' 0-^) 

If a 2 = 1, it describes an Hamiltonian wave system which is exactly solvable but not 
Liouville integrable. It exhibits a degenerative dispersion, the asymptotic states for t — > 
±oo do not coincide and an infinite number of invariants of motion exist. If a 2 = — 1, 
it describes an Hamiltonian wave system exactly solvable and completely integrable. It 
exhibits non-degenerative dispersion and lack of decay, and the asymptotic states for t — > 
±oo coincide upon imposing rapid-decrease boundary conditions. Moreover, the number 
of invariants of motion remains infinite [8] . 

The general set of KP equations may be described by using the first-order operator 

r-£-4 (1.10) 
with x G R, and the associated Lax pseudo-differential operator 

oo 

L = T + ^2u k (x,t u t 2 ,t 3 , ...)T~ k (1.11) 
k=i 

where the functions u k are here called the 'potentials'. By doing so, one allows in general 
for their dependence on an infinite number of time variables t = (t±, t 2 , ...,t p , ...). On 
assuming that T _1 is a well defined inverse operator (see section 2), so that 

T ■ T -1 = T~ l T = 1 (1.12) 

one can compose the Lax operator with itself, giving rise to its 'powers', i.e. L n = L-L n_1 , 
for all n = 2,3, oo. Each such power has a differential part, denoted by B n . To begin 
one sets 

B 1 = T (1.13) 
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and, by virtue of (1.12), one finds 

B 2 = T 2 + 2 Ul (1.14) 

B 3 = T 3 + 3 Ul T + 3(u 2 + u hx ) (1.15) 

and so on. The KP hierarchy of integrable equations is then defined by the generalized 
Lax equation [9, 10] 

dL 

— = [B n ,L] = B n L-LB n (1.16) 

and by the Zakharov-Shabat equation 

dB m dB 



dt dt 



[B n ,B m ] (1.17) 



which may be seen as the compatibility conditions of the linear equations 

Lip = Xip (1.18) 

and 

w n = B »* < L19 > 

for all n, under the assumption that 

° X = 0. (1.20) 



dt n 

At this stage, since t\ plays the same role as x, t\ or x are used without distinction in the 
literature [10]. Once the equations (1.16) are written for all values of n, the coefficients of 
T~ k are equated, and this leads to an infinite set of equations 

= Vkn (1-21) 



dt 



n 



where tpkn are certain differential polynomials in the potentials and their derivatives. For 
example, from the equations [10] 

du\ 

ST = Ul,xx + 2u2, x (1-22) 



dt 2 

dui 



du2 

= u 2 , xx + 2u 3jX + 2mu ljX (1-23) 



... =u liXXX + 3u2, xx + Su 3iX + 6u 1 u ljX (1-24) 
013 

one gets the KP equation for u\. 

In section 2 we derive the kernel of the Lax pseudo-differential operator by using a 
Green-function method and the theory of distributions when all potentials only depend 
on x. In section 3 we obtain the symbol from the kernel of section 2, and derive suitable 
majorizations which ensure ellipticity of the Lax operator. Strong ellipticity is studied in 
section 4. Behaviour of the ellipticity conditions under KP flows is investigated in section 
5, and concluding remarks are presented in section 6, while the appendix describes relevant 
details. 

2. The Lax operator and its kernel 

Following [10], we first consider a 'restricted' form of the Lax operator, for which the 
potentials Uk only depend on x. The general form (1.11) will be restored in section 5, 
where time evolution is studied (cf section 3 of [10]). 

Once the operator (1.10) is given, the inverse operator T -1 is an integral operator 
with kernel given by the Green function Gi(x, y) of T. Its action on any function / in its 
domain reads 

/oo 
G 1 (x,y)f(y)dy (2.1) 
-00 

where the Green function G\ obeys the equation 

T x G 1 (x,y) = S(x,y). (2.2) 
More precisely, the Green function Gi is a kernel which solves the equation 

■^G 1 (x,y) = V x + y (2.3) 
6 



and the jump condition [11] 

lim G x {x,y)- lim G 1 (x,y) = l. (2.4) 

x^y+ x^y~ 

The problem described by Eqs. (2.3) and (2.4) is solved by 

G!{x,y) = A hl (y) if x > y (2.5a) 

G 1 (x,y) = A 1A (y)-l if x<y (2.56) 

where Ai^(y) is an arbitrary smooth function of y unless a suitable boundary condition is 
specified (see below). 

Similarly, the operators T~ 2 , T~ 3 and so on are integral operators with kernel given 
by the Green function of T 2 , T 3 , . . . , respectively. For example, the operator T 2 = T T has 
a Green function G2(x,y) satisfying the differential equation 

° G 2 (x,y) = V x^y (2.6) 



dx 2 



the continuity condition 



and the jump condition 



lim G 2 (x,y)= lim G 2 (x,y) (2.7) 

x^y+ x^y~ 



lim ^ - lim ^ = 1. (2.8) 

x^y+ OX x^y- OX 

Equations (2.6)-(2.8) are solved by 

G 2 {x, y) = A h2 (y) + A 2 , 2 (y)x if x > y (2.9a) 

G 2 (x, y) = y + A h2 (y) + (A 2 , 2 (y) - l)x if x < y (2.96) 

where now two arbitrary functions Ai t2 and A 2t2 are involved because G 2 is the Green 
function of a second-order differential operator. 
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It is therefore clear that, assuming for the time being that Uk only depends on x, the 
Lax operator (1.11) can be viewed as an integral operator whose action is given by 



/oo 
K(x,y)^(y)dy (2.10) 
-oo 

with kernel 

oo 

K(x, y) = -5'{x, y) + ^2 u p( x ) g p (x, V) (2- 11 ) 

P =i 

where we have used the well known distributional action of the first derivative of the Dirac 
delta functional [12], and the Green function G p (x,y) can be expressed in the form 

p 

G p (x,y) = J2Cr, P (y)x r ~ 1 - (2-12) 



r=l 

The 'coefficients' C r , p are actually functions of y obeying a law of the type (see (2.5) and 
(2.9)) 

c r , P (y) = A r, P (y) if x > y (2.13a) 

C r , p (y) = B rjP (y) if x < y (2.136) 

where the coefficients B r ^ p {y) in (2.13b) can be expressed in terms of the coefficients A rjP (y) 
after imposing the continuity conditions 

d q G 8 q G 
lim + ^- lim ^=0 V , = 0,l,...,p-2 (2.14) 

x^y+ OX q x^y- OX q 

and the jump condition 

l im ^JfP - l im = 1. (2.15) 

x^y+ dxP- 1 x^y- dxP- 1 V ' 

In [7], the Green function Gi(x,y) given by (2.5a) and (2.5b) has been written in the 

form 

Gi (x, y) = l - [6(x -y)- 6(y - x)] (2.16) 



where 9 is the step function such that 9(x) = 1 if x > 0, 9(0) = \, 9(x) = if x < 0. 
Equation (2.16) corresponds to choosing 



^i.i(y) = g 



in Eqs. (2.5a) and (2.5b). The operator T 1 is then the integral operator 



f(y)dy- - / f(y)dy 



(2.17) 



(2.18) 



Similarly, the operator T 2 turns out to be the integral operator 



OO J — OO 



X 



T •-' : ./' >\ f dy /'" /!:!</; - J / dy / /( -.)■■/- 



-co Jy 

oc 



1 /-oo t-y i ,-oo ,-oo 

j dy f(z)dz+- / dy / /(z)cte. 



(2.19) 



By comparison with (2.9a) and (2.9b) this leads to the evaluation of A lj2 (y) and A 2j2 (y), 
and the procedure can be iterated (in principle) to obtain all A rtP (y) coefficients in (2.13a), 
while the B rjP (y) are obtained after imposing (2.14) and (2.15) as we said before. The 
details of the construction are indeed a bit involved, and hence it is worth showing what 
can be done with the integral operator T~ 2 given in (2.19). On the one hand, Eqs. (2.9a) 
and (2.9b) lead to 



(T- 2 f)(x) = f dy\A lt2 (y) + xA 2 , 2 (yj\f(y) 

J — oo 



+ 



dy 



y + A h2 (y) + x(A 2 , 2 (y)-l) f(y). 



(2.20) 



On the other hand, by virtue of (2.19), (T 2 f)(x) is also given by 



f x 1 f°° ( 1 

(T- 2 f)(x) = dy-h(y) + J dy (-^(y) 



(2.21) 



where 



/y poo 
f(z)dz- / f(z)dz. 
-oo J v 



(2.22) 



Direct comparison of the representations (2.20) and (2.21) of (T 2 f)(x) yields therefore 
the equations 



A i .2(!j) + -rA,o(!j) f(y) = \h{y) 



y + A lj2 (y) + x(A 2j2 (y) - l)j f(y) = --%). 



(2.23) 
(2.24) 



The addition of (2.23) and (2.24) leads to 



2A 1>2 (y) + 2xA 2j2 (y) +y-x f(y) = 



(2.25) 



which is satisfied for all f(y) if and only if 



Such a system is solved by 



2A lj2 (y) + y = 



(2A 2)2 (y) - l)x = 0. 



MM = 



1 



(2.26) 
(2.27) 

(2.28) 
(2.29) 



which provides the desired explicit formula for the Green function G 2 (x, y), upon insertion 
into (2.9a) and (2.9b). 

3. Symbol and ellipticity 

Recall now that, if L is a pseudo-differential operator defined by a kernel K, this is related 
to the symbol p(x, £) by the equation [3] 



K(x,y) = (2 



n)- ( 



j(x-y)-£ 



(3.1) 



This equation can be inverted to give a very useful formula for the symbol, i.e. (cf Eq. 
(2.1.36) in [3]) 

p(x,£)= I e~ lz< K{x,x- z)dz. (3.2) 
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Equation (3.2) is a key formula for our investigation, because the ellipticity of L is defined 
in terms of its symbol, as we know from the introduction, following [1]. 
In our problem, which involves x G R, the integral (3.2) reduces to 



p(x 



/oo 
e- iz ^K(x,x- z)dz 
-oo 



(3.3) 



where the kernel K{x,y) is expressed by (2.11)-(2.13), and we have to check that the 
inequality (1.7) is satisfied for |£| > Co to obtain ellipticity. Indeed, the symbol (3.3) turns 
out to be 



P{x,0 = 



e-^S f {z)dz 



/oo oo P 

e -iz£, ^2 u p (x) ^2 C r , P (x - z)x r ~ 1 dz 
-°° p=l r=l 



(3.4) 



and hence obeys the inequality 

\p(x,0\> 



*8'{z)dz 



/oo oo P 

e -iz£ U p (x) Cr,p( X ~ Z)x r ~ 1 dz 
-°° p=l r=l 



(3.5) 



Of course, the first integral on the right-hand side of (3.5) becomes meaningful within the 
framework of Fourier transform of distributions [13]. In the simplest possible terms, one 
has actually to consider the parameter-dependent integral (here a > 0) 



AO = J 



e -az e -iztS'(z)dz 



/oo 
S(z)(-2az - ii)e~ az2 - l ^dz. 
-oo 



(3.6) 



By virtue of the property defining the Dirac delta functional, according to which [12] 



(SJ) = f(0) 



(3.7) 
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the integral (3.6) equals — z£, and hence the first term on the right-hand side of (3.5) equals 
|£|. Now we distinguish two cases, depending on whether 



f{*,t) = \t\ 



is positive or negative. If 



e izj u p (x) C r , p (x — z)x r 1 dz 
P =i 



r=l 



(3.8) 



/(*,£)> 



(3.9) 



holds, the majorization (1.7) for the ellipticity of the restricted Lax operator is satisfied 
provided that, for |£| > C , 



/oo oo p 

e - l2 £ ^ u p (x) ^ C r>p (x - z)x r ~ 1 dz 
-°° p=l r=l 



^lei-cr^i + KD' 



(3.10) 



In contrast, if 



f{x,t)<0 



(3.11) 



holds, the restricted Lax operator is elliptic provided that 

p 



/OO — ' f 
e -iz£ ^2 u p (x) ^2 C r , P (x - z)x r ~ 1 dz 
-oo i i 



p=l r—1 

>c -cr\i+m 



>\z\-c^(i + m 



(3.12) 



for |^| > Co. If the order d of the Lax operator is positive, we can further write that, for 
|C| > Co, the majorization (3.10) becomes 



e izj y u p (x) ^ C r ^ p {x — z)x r x dz 
P =i 



r=l 



< iei - c^cg 



(3.13) 



If we are interested in sufficient conditions we can point out that, since the inequality 
(3.5) is always satisfied, whereas (1.7) only holds when L is elliptic, the sufficient condition 
for ellipticity of the restricted Lax operator is expressed by 



cr 1 (i + \c\) d <\c\ 



/oo CO p 

e -tz$ ^2 Up (x) ^2 C r , P (x - z)x r ~ 1 dz 
-°° p=l r=l 



(3.14) 
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for |£| > Co. This leads in turn to the inequality (3.10). 

To sum up, if the function / : T*(R) — > R defined by (3.8) has no zeros (which 
is already a non-trivial requirement on the potentials u p ), the restricted Lax operator 
is elliptic provided that either (3.10) or (3.12) is satisfied. The majorization (3.10) is 
further simplified in the form (3.13) in case of positive order of the Lax operator. A 
sufficient condition for ellipticity is given instead by (3.14), which coincides with (3.10). 
In particular, when |£| = Co and the equality sign is chosen in (3.14), the order d of the 
restricted Lax operator can be evaluated by the formula 



d = 



log(Ci(C -Jg)) 
log(l + C ) 



(3.15) 



where 



/oo oo P 

e -iz£ ^2 Up (x) ^2 C r , P (x - z)x r ~ 1 dz 
-°° p=l r=l 



(3.16) 



h = SU Px 6 R 

bearing in mind that the modulus of £ equals the constant Co- 
4. Strong ellipticity 

In a thorough analysis of the ellipticity properties, strong ellipticity should also be studied. 
For this purpose, following [3] , we assume that the symbol of the restricted Lax operator 
is polyhomogeneous, in that it admits an asymptotic expansion of the form 



(4.1) 



1=0 



where each term p^-i has the homogeneity property 



Pd-l{x,l£,) = l d l p d -i(x,0 



(4.2) 



for t > 1 and |£| > 1. The principal symbol p° of the Lax operator is then, by definition, 



p°(x,£) =p d (x,£). 



(4.3) 
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Strong ellipticity is formulated in terms of the principal symbol, because it requires that 



Rep (x,O = l[p (x,O+P°(x,O*] >c(xm d (4.4) 

where x G R, c(x) > and |£| > 1. In other words, given a positive function c, the 
product c(x)|£| d should be always majorized by the real part of the principal symbol of 
the restricted Lax operator. Indeed, the symbol (3.4) is such that 

p(x,7f) = 

/•oo oo p / \ 

+ 7" 1 / e"** ^ «p(x) ^ C r , P ( x - - x r ~ x dz. (4.5) 

p=l r=l ^ 7 ^ 

By virtue of (4.1), (4.2) and (4.5) we find that 

/■oo oo v , n 

- i 7 £ + 7" 1 / e- te « ]T tip(x) Z" 1 ^ 

"'- 00 p=l r=l ^ 7 ^ 

CXD 

Moreover, the term on the right-hand side of (4.6) with / = should be the one occurring in 
the condition (4.4) for strong ellipticity. A mathematical advantage of strong ellipticity lies 
in the possibility of having a well defined functional trace of the heat semigroup associated 
to the Lax operator [1, 3]. 

5. Behaviour of the ellipticity conditions under KP flows 

To study the preservation (or violation) of the ellipticity conditions under KP flows one 
has to analyze the following problem: suppose that the conditions (3.10) or (3.12) are 
satisfied for t = 0. Are they still valid for all or some t > 0? 

This means that we consider again the potentials Uk as in Eq. (1.11), i.e. as functions 
of x and t, where t is a concise notation for infinitely many time parameters (t\,t2, ■■■). 
It should be stressed that we consider only one spatial variable and infinitely many time 
parameters, since otherwise it would be problematic, at least for the authors, to generalize 
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formulae like (3.2) aimed at obtaining the symbol from the kernel of the operator. For our 
purposes it is convenient to use formulae generating such potentials by means of a single 
function. This is made possible by the r- function (see appendix), because one finds [10] 



d 2 



1 / d s d 2 
t) = -[ -—^ + — — ) logr (5.2) 



2 V dx 3 dxdt 2 
1 / <9 4 d 3 d 2 \ 

and infinitely many other equations of the general form 

u k (x;t) = F fc (logr) (5.4) 

where F k is, in general, a non-linear function of logr. Equations (5.4), with k ranging 
from 1 through oo, should be inserted into the ellipticity conditions (3.10) and (3.14), 
substituting therein u k with ^(logr) for all k. The resulting majorizations involve non- 
linear functions of the logarithm of the r-function. 

Further progress can be made by considering a 'truncated' Lax operator, e.g. 

L = T + u 1 (x;t)T~ 1 + u 2 (x;t)T- 2 . (5.5) 

This should not seem an arbitrary simplification, because the Lax operator is obtained 
from the W operator of the appendix as [10] 

L = WTW-\ (5.6) 

Now both W and its 'truncated version' [10] 



Wm = l + J2^kT- k (5.7) 



fc=i 
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satisfy the Sato equation (A8), from which the generalized Lax equation (1.16) is eventually 
obtained. Thus, operators like L in (5.5) can be obtained from (5.6) if W is replaced by 
W m therein. In this case, on defining 

2 p 

F(x, z; t) = u p (x; t) C r}P (x — z)x r ~ 1 
P =i 



ui(x; t)C ltl (x — z) + u 2 (x; t) Ci )2 (x - z) + xC 2)2 (a: - 2;) 



(5.8) 



the integral on the right-hand side of the ellipticity condition (3.14) reduces to (Ji + J 2 + 
J 3 )(£,£;t) where, bearing in mind that (see (2.17), (2.28) and (2.29)) 



one finds 



1 1 

Ci,i(x -z) = C 2 , 2 (x -z) = -i£z>0, -- if z < 



C 1>2 (x-z) = -| + |ifz>0, I -|ifz<0 



/oo 
e~ l2 ^wi(x; t)Ci ) i(x - 
-00 



= iti uc;t) hm -e 2 sin — 
v y 6^00 2 



(5.9) 
(5.10) 



(5.11) 



/oo 
e" i2C « 2 (x; t)C 1:2 (x - z)dz 
-00 

/ .^1 . / , v d Ji(x,£;t) 

= — xm 2 (x;cJ r — \-U2\x;t)i- 



ui(x;t) 



<i£ ui(x;t) 



(5.12) 



/oo 
e~ lz? 'U 2 (a;; t)xC 2)2 (x - 
-00 

Ji(x,£;*) 



XU2(X] t)- 



ui(x;t) 



(5.13) 



In these equations, the infinite upper limit of integration can be recovered by taking the 
limit as b — > 00 of integrals from to 6 (the lower limit being amenable to by virtue of 
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(5.9) and (5.10). However, divergences remain, and hence we are only able to obtain well 
defined formulae by integrating up to finite values of b. The cancellation of terms involving 
xu 2 (x; t) is thus found to occur on performing the sum, and our integral reads eventually 

( Ji + J 2 + J 3 ) b (x, £; t) = Ul (x; t)F b (£) + u 2 (x; t)i-^F b (£) (5.14) 



having defined 



F b (0 = ^e- ib hm b -^. (5.15) 



Now we point out that 

\ Ul (x;t)F b (Z) + u 2 (x;t)iF{,(Z)\ < \ Ul (x; t)F b (£)\ + \u 2 (x; t)iF^)\ 



<|M*; t )| + \u 2 (x;t)\. (5.16) 



Thus, a sufficient condition for the validity of the majorization (3.14) is expressed by 

ICI > Cr\l + \Z\) d + ^Mx;t)\ + Q + |) \Mx;t)\- (5.17) 

It should be stressed that divergent integrals occur already in the 'time-independent' el- 
lipticity condition (3.14). When the time parameters are introduced, it may be easier or 
harder to fulfill ellipticity, depending on the behaviour of ui,u 2 ,... (which, in turn, all 
depend on the r- function). 

6. Concluding remarks 

Our paper has been motivated by the need to obtain a deeper understanding of the basic 
structures of modern non-linear physics. It is indeed well known that the Sato equation 
(A8) generates the generalized Lax equation (1.16), the Zakharov-Shabat equation (1.17) 
and the inverse spectral transform scheme [10]. Moreover, an infinite number of nonlinear 
evolution equations (i.e. the KP hierarchy), of which the KP equation is the simplest 
nontrivial one, share solutions, and the r-function makes it possible to express all such 
solutions. 
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The work in the present paper is the first step towards a rigorous investigation of 
the ellipticity properties of the Lax pseudo-differential operator. We have found that, to 
achieve ellipticity, including its strong form, the various potentials Uk(x) are no longer 
arbitrary, but should be chosen in such a way that the following conditions hold: 

(i) The function / : T*R — > R defined in (3.8) has no zeros. 

(ii) The majorization (3.14) holds (some care is actually necessary to deal with the integrals 
on the right-hand side of (3.14), as is clear from the analysis performed in section 5). 

(iii) The asymptotic expansion (4.6) can be obtained. Note, however, that violation of 
(4.2) (i.e. lack of homogeneity) for |£| < 1 can cause logarithmic terms in the asymptotic 
expansion of the kernel defined by (2.11)-(2.13) (cf [14]). 

Moreover, on allowing for the time evolution of the potentials in the Lax operator, 
now viewed as functions of x and of infinitely many time variables, we have obtained an 
explicit ellipticity condition in terms of the r-function, when attention is restricted to the 
'truncated' Lax operator (5.5). It now remains to be seen how to deal with the infinite 
sum over p in the ellipticity condition (3.14) when the potentials u p (x;t) appropriate for 
the 'full' Lax operator (1.11) are instead considered. The form of u\ and u 2 remains the 
one given in (5.1) and (5.2), but the occurrence of an infinite number of such potentials 
makes it hard to re-express the time-dependent form of the majorization (3.14). The work 
in [15] has indeed obtained a very useful formula for the r-function, but this remains of 
little help when the infinite sum over all potentials is performed. 

A further interesting issue is the investigation of ellipticity for the formulation of 
KP hierarchy considered in [16], where the main new technique, when compared to the 
traditional approach to the generalized Lax equation, consists of replacing the Lax operator 
by an n th -order formal pseudo-differential operator 

n-2 

L n = T n + 1j TJ n ^ 2 - 

j=-oo 
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The authors of [16] have been able to factorize L n into n — 1 first-order formal differential 
operators Ak, 1 < k < n — 1, and one first-order formal pseudo-differential operator A n , 
i.e. 

L n = A n A n _ 1 ...A 2 A l (6.2) 

where 

A k =T + 7] kjX l<k<n (6.3) 

n 

X>M = (6.4) 

fc=i 

-l 

A n = A n + Ki T '- (6-5) 

j = -oo 

The results and unsolved problems described so far seem to show that new exciting 
developments might be obtained from the effort of combining some key techniques of non- 
linear physics and the tools of linear and pseudo-differential operator theory. In particular, 
the mathematical requirement of ellipticity in the various forms considered in sections 3- 
5 restricts the potentials Mfc in a form not previously considered in the literature to our 
knowledge, which might be used to select the realizations of the Lax operator one is 
interested in. 
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Appendix 

Since the general reader is not necessarily familiar with the theory of r-functions, we 
summarize the main properties hereafter. Following [10], we study for the operator W m 
defined in (5.7) the ordinary differential equation 

W m d m f{x) = (d m + w l (x)d m ~ 1 + ... + w m (x))f(x) = (Al) 
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which has m linearly independent solutions f^\x), f( m \x). On writing equation (Al) 
m times with / = f^\x), / = /( m )(x), one finds a linear system which can be solved 
for wj(x), for all j = 1, ...,m, and hence W m is found from the definition (5.7). 

When the Wj are taken to depend also on infinitely many time variables (t±, t p , ...), 
the operator W m (x; t) is found to be 



W m {x;t) = 



detA 
r(x;t) 



(A2) 



where, given the functions Kq (x; t) and hn'(x;t) such that 



d d r 



dt n dx n 



h { j) (x;t) = n=l,2, 



(A3) 



hg\x;0) = fW(x) 



hU\x;t) 



one has (see the definition (1.10)) 



9 h^(x;t) = ^\x;t) 



dt n 



and 



(A4) 
(A5) 







ILq 


r-p—m ^ 




A EE 




ui.ru) 
••• a m-l 
u(m) 


T -l 

i y 


(A6) 


r(x;t) 


ee det | 


'tf> ... 
\ U m-l 




(A7) 



The function r is the r- function used in Eqs. (5.1)-(5.4), and the time evolution of the 
operator W m (x;t) is determined by the Sato equation 



dt n 



"' B n W m - \V,„T" 



(A8) 



where the operators B n are the same occurring in the generalized Lax equation (1.16). 
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